C/3 

x. 



Green's functions for Neumann boundary conditions 



Jerrold Franklii 

Department of Physics, Temple University, Philadelphia, PA 19122-6082 

(Dated: August 28, 2012) 

Abstract 

Green's functions for Neumann boundary conditions have been considered in Math Physics and 
Electromagnetism textbooks, but special constraints and other properties required for Neumann 
boundary conditions have generally not been noticed or treated correctly. In this paper, we derive 
an appropriate Neumann Green's function with these constraints and properties incorporated. 
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I. INTRODUCTION 



The Green's function method for solving Sturm-Liouville problems of the form 

£ r 0(r) = V ■ [p(r) V0(r)] + g(r)0(r) = p(r) (1) 

is described in many textbooks. We list some of the more recent texts that treat Green's 
functions in references 1-4, and review the usual textbook treatment below. For simplicity, 
we will take q(r) = 0. 

A motivation for using a Green's function, G(r, r'), is that it satisfies homogeneous bound- 
ary conditions that makes it easier to solve for than the original problem with inhomogeneous 
boundary conditions. Application of Green's theorem 

J [0(r')£ r ^(r, r') - G(r, r / )£ r ^(r / )]dr' = f dS' ■ W)p{v') V'G(r, r') - G(r, r')p(r') V'0(r')] 

(2) 

provides the solution to the original problem if G(r, r') satisfies 

£ r ,G(r,r') = 5(r-r'). (3) 

In that case, Eq (j2J) leads to 

0(r) = | G{r, x')p{r')dr' + J dS' • [0(r>(r')V'G(r, r') - G(r, r>(r')V'0(r')]. (4) 

The two simplest boundary conditions for which the Green's function method is applicable 
are the Dirichlet boundary condition for which the solution 0(r) is given on all bounding 
surfaces, and the Neumann boundary condition for which its normal derivative n • V0(r) is 
given. The Dirichlet Green's function is generally used for electrostatic problems where the 
potential is specified on bounding surfaces, while the Neumann Green's function is useful for 
finding temperature distributions where the bounding surfaces are heat insulated or have 
specified heat currents. 

For the Dirichlet boundary condition, the Green's function satisfies the homogeneous 
boundary condition 

Go(r, r') = 0, for r' on all bounding surfaces. (5) 
This reduces Eq. (TjJ to 

0(r) = J G D (v, r')p{y')\dr' + J dS' • W)p{v')V'G D {v, r% (6) 
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which is the solution to the Dirichlet problem once Gi)(r, r') is known. The Dirichlet Green's 
function is the solution to Eq. ([3]) satisfying the homogeneous boundary condition in Eq. 

©■ 

Textbooks generally treat the Dirichlet above, but do much less with the Green's 

function for the Neumann boundary condition, and what is said about the Neumann case of- 
ten has mistakes of omission and commission. First of all, the Neumann boundary condition 
for the solution <p(r) must satisfy the constraint 

|dS.p(r)W(r) = Jpdr, (7) 

which follows from applying the divergence theorem to Eq. ([1]). Most texts do not mention 
this important constraint on the Neumann boundary condition. 5 

There are cases where the boundary condition is Neumann on some surfaces and Dirichlet 
on others. In those cases, the normal derivative of <fi on the Dirichlet surfaces automati- 
cally adjusts to satisfy the constraint surface integral. But for pure Neumann boundary 
conditions, the normal derivative must satisfy the constraint or no solution exists. In the 
following sections, we will assume pure Neumann boundary conditions for which constraint 
equation (J7j) holds. We treat the one dimensional Neumann Green's function in Section 2, 
and then the three dimensional case in Section 3. 



II. ID NEUMANN GREEN'S FUNCTION 

We first review how the Dirac delta function arises when a function f(x), defined in the 
finite range < x < L, is expanded in orthonormal Dirichlet eigenfunctions u n (x) of a 
Sturm-Liouville operator C x . The eigenfunctions satisfy the differential equation 

= \ n u n {x). (8) 

with homogeneous Dirichlet boundary conditions 

u n (0) = 0, u n {L) = 
The function f(x) can be expanded in the Dirichlet eigenfunctions as 

oo 

f( X ) = b nUn(x), (10) 
n=l 
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dx 



p(x 



du n 


dx 



(9) 



with the expansion coefficients b n given by 

u* n (x)f(x)dx. (11) 



o 

If Eq. f lTTj) for the expansion coefficients is substituted into Eq. (fTUj) for f(x), and the sum 
executed before the integral, we get 

f(x)= / dx'Y,<(*>n(x)f(x'). (12) 
J ° n=l 

From the definition of the Dirac delta function by its sifting property, 

f(x) = [ L dx'6(x-x')f(x') (13) 

for x and x' in the range [0, L], we see that the delta function can be represented by a sum 
over Dirichlet eigenf unctions as 

00 

5(x — x) = u* n (x')u n (x). (14) 



n=l 



For example, for the simple case 



£ x M n = u' n = X n u n , u n = ^2/Lsm(nnx/L), X n = - ( — J , (15) 
the delta function is represented by 

5(x — x') = {2 /Li) sm(n7rx' /L) sm(mrx/L) (16) 

n=l 

A Dirichlet Green's function that satisfies the differential equation 

C x 'Gd(x,x') = 5(x — x'), (17) 

and satisfies the homogeneous Dirichlet boundary conditions in the variable x' can be formed 
from the Dirichlet eigenf unctions as 

G D (x,x') = jr <{x ' )Un{x) . (18) 

n=l 

Acting on this Green's function with the Sturm Liouville operator £ x > removes the denomi- 
nator in Eq. ( I18p leaving the delta function of Eq. ( 1T7|) . showing that this is the appropriate 
Dirichlet Green's function. 

The above straightforward derivation for Dirichlet boundary conditions is given in most 
texts, but a corresponding derivation for Neumann boundary conditions is generally absent. 



Neumann eigenfunctions, v n (x), satisfy the same differential equation (jSJ) as the Dirichlet 
eigenfunctions, but have the boundary conditions 

<(0) = 0, v' n (L) = 0. (19) 

The expansion in Neumann eigenfunctions has a constant eigenfunction corresponding to 
a zero eigenvalue, so the expansion is given by 

oo 

f(x) = a + a n v n (x), (20) 

n=l 

The Neumann expansion coefficients a n are given by the integrals 

a n = [ 
Jo 

and 



v^(x)f(x)dx, n > 1, (21) 



1 f L 
ao= L 



[ L f(x)dx=<f>, (22) 
Jo 



where < f > represents the average value of the function f(x) over the interval [0, L\. Now 
putting Eqs. (|2"T1) and (1221) into the expansion Eq. (1201) results in 



f(x)= / dx 







oo 

^ n=l 



/CO, (23) 



so the representation of the delta function in terms of Neumann eigenfunctions is 

-r oo 

S(x-x') = T + ^v* n (x')v n (x). (24) 

L n=l 

The additional constant term 1/L is not generally recognized in textbooks. For the simple 
case of £ x v n = t>" = A n f with the eigenfunctions satisfying the homogeneous Neumann 
boundary conditions of Eq. ( TT9|) the delta function is represented by 

1 2 00 

5(x — a/) = — + — ^2 cos (nn x / L) cos(n7rx / L) . (25) 

L L n =l 

A Neumann's Green function can be formed using Neumann eigenfunctions of the oper- 
ator C x as the sum 

GnM = ^<M^M. (26) 

n=l 

This Green's function satisfies the homogeneous Neumann boundary conditions 

d x ,G N (x,x')\( x , =0 ) = 0, dz/G/v^aOlo^L) = 0. (27) 



However, because of the constant term 1/L in Eq. (I21I) . the operation on Gn(x,x') by C x > 
is 

^/Gjv^x') = - 1/L. (28) 

Thus, the Neumann Green's function satisfies a different differential equation than the 
Dirichlet Green's function. 

We now use the Green's function G^(x, x') to find the solution of the differential equation 

df 



CJ{x) = T X 



V{x 'dx 



p(x), (29) 

with the inhomogeneous Neumann boundary conditions 

f'(0) = f Q , f(L) = f L . (30) 

The boundary values must satisfy the constraint 

p(L)f' L -p(0)ti= f L p(x)dx, (31) 
Jo 

which follows from a first integral of Eq. (T29]) . The Neumann Green's function must also 



satisfy this constraint, which it does because the right hand side of Eq. (128]) integrates to 
zero. 

Green's theorem in one dimension (or integration by parts) for this differential equation 
leads to 

f L [f(x')C x/ G N (x,x') - G N (x,x')C x/ f(x')}dx' = -G N (x,L)p(L)f L + G N (x,0)p(0)f^. (32) 
Jo 

We have used the boundary conditions (T2"T|) to eliminate terms containing d x 'G(x, x') at the 
endpoints. Then, using Eq. (1281) . we get 



f(x) =<f> + [ L G N (x, x')p(x')dx' - G N (x, L)p(L)f' L + G N (x, 0)p(0)^, (33) 
Jo 

which constitutes the solution to the Sturm-Liouville problem for Neumann boundary con- 
ditions. The constant < / >, the average value of f(x), arises when the term 1/L in Eq. 
(Hi is substituted into Eq. (152]) . 

Actually, any constant can be added to the solution f(x) since the solution of the Neu- 
mann problem is only unique up to an additive constant. The constant term in Eq. (1331) 
will always be the average value of the solution because the variable terms have zero average 
value. 



Although we have used an expansion in eigenf unctions to give an heuristic derivation of 
the Green's function for Neumann boundary conditions, any function satisfying the defining 
equation ( 1251) with the boundary conditions of Eq. (127)) will be a suitable Green's function 
Gn(x.x'). For example, a Green's function for the problem 



/"(*)=*, /'(0) = 4 /' (L) = f' L = ft + L 2 /2, (34) 



is given by 



x' 2 

< x' < X : Gpfi(x, x') = — jry + X 

x' 2 

x < x' < L : G N2 (x, x') = -— + x'. (35) 

The term —x' 2 /2L in Gm and Gjv 2 provides the —1/L term in d 2 ,Gx(x,x'), and also 
satisfies the homogeneous Neumann boundary condition at x' — 0. The a;' term in Gjv2 
satisfies the homogeneous Neumann boundary condition at x' = L, and also provides a unit 
step in Gn at x' — x so that the next derivative will give the delta function in d 2 ,GN{x,x'). 
Finally, the term x in Gni makes Gn{x, x') continuous at x' = x. Putting this Green's 
function into Eq. (133)) gives the solution to Eq. (1341) . 

We note that the Green's function in Eq. (1351) is not symmetric with respect to x and 
x' . To show that a Neumann Green's function need not be symmetric, we repeat the usual 
proof of symmetry here. Applying Green's theorem to two Green's functions, G(xi, x') and 
G(x2,x'), of a Sturm-Liouville operator C gives 



\ [G(xi, x')C x iG(x2, x') — G(x2, x')C x >G(xi, x')\ 
Jo 



= [G(x 1 ,x')p(x')d x >G(x 2 , x') - G(x 2 , x')p{x')d x ,G{x\,x')] x xl J^ (36) 

The right hand side of Eq. ( 136)) vanishes for either Dirichlet or Neumann homogeneous 
boundary conditions. For a Dirichlet Green's function, with C x 'Gd(x, x') = 5(x — x'), Eq. 
(136)) reduces to 

G D {xi,x 2 ) - G D (x 2 ,x 1 ) = 0, (37) 

so a Dirichlet Green's function must be symmetric. However, for a Neumann Green's func- 
tion, C x 'Gn(x, x 1 ) = —1/L + 5(x — x'), and Eq. ( 136)) reduces to 

G N (xi,x 2 ) - G N (x 2 ,x 1 ) = — I [Gn(xi,x) - G N (x 2 ,x')]dx' } (38) 
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so a Neumann Green's function is not required to be symmetric. 

We do know however, from the eigenfunction expansion in Eq. (126]) that any real Neumann 
Green's function can be made symmetric. For instance, adding the term —x 2 /2L to the 
Green's function in Eq. f[3"5"]) will make that Neumann Green's function symmetric without 
changing any of its actions. There is no need to do this however, since the non-symmetric 
Green's function is simpler, and either form will solve the original differential equation. 

III. 3D NEUMANN GREEN FUNCTION 

In this section, we extend the one dimensional results of the previous section to three 
dimensions. In three dimensions, we seek the solution of the differential equation 

£ r 0(r) = V ■ [p(r)V0(r)] = p(r) (39) 

with the inhomogeneous Neumann boundary conditions on any bounding surface 

n-V<Kr) = /(r 5 ), (40) 

where n is the outward normal vector to the surface and f{rs) is an almost arbitrary function 
specified on all surfaces. A solution to Eq. ( 139]) exists only if the boundary conditions satisfy 
the constraint 

I dS • p(r)V0(r) = J p{r)dr. (41) 

This constraint follows by applying the divergence theorem to Eq. (I3"9l . 

The Neumann Green's function for this problem satisfies the differential equation 

£ r ,G N (r,r') = 5(r-r')-l/V, (42) 

with the homogeneous boundary condition 

n'- VG N (r,r') = (43) 

on all surfaces. This Green's function automatically satisfies the constraint 

JdS'-p(r')VG N (r,r') = 0, (44) 

which follows by applying the divergence theorem to Eq. fl42]) . 
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The extra term 1/V in Eq. (1421) is the 3D equivalent of the term 1/L in Eq. (|28|) . and 
arises due to the constant term in any expansion using Neumann eigenfunctions. The 1/V 
goes to zero for an infinite volume (the so called 'exterior problem'), but is important if the 
volume considered is finite (the 'interior problem'). 

The solution to Eq. (|3"9"|) is given by 

0(r) =< > + J G(r, r')p(r')dr' - J dS' • G(r, r')p(r') V'0(r'), (45) 

which follows from Green's theorem, and Eqs. (j42p and (I4"3"j) for GAr(r, r'). As in the ID 
case, < > is an arbitrary constant that equals the average value of <fi in the volume. It 
can be chosen to be zero to simplify the equation. 

An alternate method to develop a Neumann Green's function has been proposed in Ref. 4. 
That method keeps Eq. ([3]) for the action of C r r on the Greens function. The normal 
derivative of the Green's function must then satisfy the constraint 

p(r')n- V'G(r,r') = 1/S (46) 

on all bounding surfaces, where the constant S is the total area of the bounding surfaces. 
This constraint follows from applying the divergence theorem to Eq. ([3]). 

This method does not always work. Kim and Jackson^ have applied it to the case of 
two concentric spheres, where spherical symmetry allows the constant normal derivative in 
Eq. f|46l) . while also satisfying Eq. ([3]). However, for more general geometries, such as for the 
temperature distribution inside a rectangular parallelepiped with insulated walls, there is 
no function that satisfies Eq. ( 146|) . while also satisfying Eq. ([3]). To satisfy Eq. ( )46l) at both 
the left and right hand face of the parallelepiped would require either a periodic function 
or a function with a discontinuous derivative, neither of which would be compatible with 
Eq. 03]). On the other hand, Eq. f )26|) provides a suitable Neumann Green's function for the 
parallelepiped, or for any geometry that allows a complete set of Neumann eigenfunctions. 
This alternate method would also not work for any one dimensional problem, because it 
would require the first and second derivatives of G(r,r') to each be of dimension 1/length, 
which is not possible. 
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IV. SUMMARY 



To summarize, we have found seven essential differences between the Neumann Green's 
function treated here and the Dirichlet Green's function which is generally treated in Math 
Physics or Electromagnetism texts for the solution of a partial differential equation of the 
form 

£ r 0(r) = V ■ [p(r)V(r)]0(r) = p(r). (47) 
We list again the different Neumann properties below: 

(1) The Neumann boundary condition for the solution <p(r) is 

n-V0(r)=/(r s ), (48) 

(2) with the constraint 

J dS ■ p(r)V0(r) = j p(r)dr (49) 
required for any solution to exist. 

(3) The Neumann Green's function satisfies the differential equation 

C rl G N (r,r') = 5(r-r')-l/V, (50) 

(4) with the homogeneous boundary condition 

n' • VG N (r, r') = (51) 

on all bounding surfaces. 

(5) The solution to Eq. ( 1471) is given by 

0(r) =< > + J G{r, r')p(r')dr' - J dS' ■ G(r, r>(x')V'0(r'), (52) 

(6) where < <fi > is an arbitrary constant that equals the average value of 0(r) in the volume. 

(7) The Neumann Green's function is not necessarily symmetric, but can always be 
made symmetric by adding a function of r to G^(r, r'). Adding any function of r to a 
Neumann Green's function does not change its actions. 
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Of these seven differences, only number (1) is generally mentioned in Math Physics or 
EM textbooks. Some books^ - — give the homogeneous boundary condition in Eq. (j5ip . but 
don't mention that it is inconsistent with omitting the 1/V term in Eq. (1501) . Reference 4 
recognizes this by making the Green's function boundary condition inhomogeneous, but this 
is not always a satisfactory remedy, as we have discussed above. 
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